Limits on the spin-dependent WIMP-nucleon cross-sections 
from the first science run of the ZEPLIN-III experiment 
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We present new experimental constraints on the WIMP-nucleon spin-dependent elastic cross- 
sections using data from the first science run of ZEPLIN-III, a two-phase xenon experiment searching 
for galactic dark matter WIMPs based at the Boulby mine. Analysis of ~450 kg-days fiducial 
exposure revealed a most likely signal of zero events, leading to a 90%-confidence upper limit on 
the pure WIMP-neutron cross- section of cr„ = 1.8 x 10"^ pb at 55 GeV/c^ WIMP mass. Recent 



calculations of the nuclear spin structure based on the Bonn CD nucleon-nucleon potential were 
used for the odd-neutron isotopes ^^^Xe and ^"^^Xe. These indicate that the sensitivity of xenon 
targets to the spin-dependent WIMP-proton interaction is much lower than implied by previous 
calculations, whereas the WIMP-neutron sensitivity is impaired only by a factor of ^^2. 
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ZEPLIN-HI has recently completed its first run at 
the Boulby Underground Laboratory (UK) searching for 
weakly interacting massive particles (WIMPs), which 
have been proposed to explain the non-baryonic cold 
dark matter in the Universe. In several supersymme- 
try (SUSY) extensions to the Standard Model of parti- 
cle physics, the lightest supersymmetric particle (LSP) is 
stable and its relic abundance could account for the mat- 
ter content observed today [1]. In most SUSY scenarios 
this particle is the neutralino, x, which is a linear com- 
bination of the SUSY partners of the electroweak gauge 
bosons (gauginos) and Higgs bosons (higgsinos). Neu- 
tralinos are expected to scatter elastically from ordinary 
matter, producing low-energy nuclear recoils, predomi- 
nantly through scalar (or spin- independent, SI) interac- 
tions. This is especially so for heavier elements (yl>30, 
with A the number of nucleons) , for which the scattering 
involves the entire nucleus rather than individual nucle- 
ons. This coherence enhancement of the scalar term in 
the interaction cross-section is proportional to A^. Our 
SI result excludes a WIMP-nucleon cross-section above 
7.7x 10"^ pb for 55 GeV/c^ mass with 90% confidence [1]. 
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Nonetheless, axial-vector (or spin-dependent, SD) in- 
teractions could dominate in some SUSY scenarios in re- 
gions where the SI term is suppressed [3i, 1^] ■ In addition, 
interpretation of the DAMA annual modulation Q in 
terms of nuclear recoils caused by a dominant SD interac- 
tion [7] has remained viable for WIMP masses below ~15 
GeV until recently. It is therefore important to continue 
pursuing the search in the SD channel. Xenon targets 
have good sensitivity to the WIMP-neutron interaction 
since approximately half of natural abundance consists 
of the odd-neutron isotopes ^^^Xe and ^^^Xe (the proton 
interaction is highly suppressed for these isotopes, and 
neither is likely in the even-even nuclei). 

It is also possible that a SD inelastic interaction with 
these nuclei could be significant for relatively heavy 
WIMPs [1, Deexcitation of low-lying nuclear states 
in i29Xe and ^^^Xe (emitting 39.6 keV and 80.2 keV 7- 
rays, respectively) provides a detection mechanism with 
low effective energy threshold, since the 7-rays can be 
efficiently detected. Although neutrons also scatter in- 
elastically producing the same signatures, differences in 
spatial distribution between signal and background might 
be exploited in very large detectors - assuming that the 
dominant 7-ray backgrounds could be suitably mitigated. 
Most direct search experiments are still focused on de- 
tecting nuclear recoils down to ^keV energies; this is so 
in the case described in this Letter. 
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WIMP search experiments propose to measure the to- 
tal WIMP-nucleus elastic cross-section, a a- Assuming 
dominance of either the SI or the SD term, one can write: 



(J A = ^G\'^\Ca 



(1) 



where Gi? is the Fermi weak-coupling constant, ^la is 
the reduced mass of the WIMP-nucleus system, and Ca 
is an enhancement factor which depends on the type of 
interaction and, possibly, the WIMP composition. For 
the SI case, Ca oc (i.e. Ca is model-independent for 
Majorana WIMPs). For a dominant SD case, Ca involves 
the total nuclear spin J instead: 

C^ = -(ap(5p)+a„(5„))'^, (2) 

TT J 



where (5'p,„) are the expectation values for the proton 
or neutron spins averaged over the nucleus and 
the effective WIMP-nucleon coupling constants. The lat- 
ter depend on the WIMP particle content (e.g. more 
higgsino- or gaugino-like) as well as the quark spin dis- 
tributions inside the nucleons. 

Two spin distribution calculations were adopted, both 
based on the nuclear shell model but using different 
nucleon-nucleon (nn) potentials (see, e.g., reviews 

Their ability to reproduce experimental measure- 
ments of the magnetic moment is the standard bench- 
mark, as this observable is reasonably similar to the 
WIMP-nucleus scattering matrix element. For compar- 
ison with other Xe experiments we use spin structures 
with the Bonn A potential for both isotopes [l^l ■ Agree- 
ment of the magnetic moment (using effective ^-factors) 
is quite reasonable: within 19% and 8% for ^^^Xe and 
"'^^^Xe, respectively. In addition, new calculations based 
on the Bonn-CD G-matrix have become available which 
improve on these figures significantly (to 1% and 2%) for 
both isotopes The spin expectation values {Sp^„) are 
given in Table I for both cases. 

In order to compare ([T]) with an experimental result 
a nuclear form factor, F'^{q), must be assumed at mo- 
mentum transfer g > to account for finite nuclear size. 
Although cross-sections are conventionally given in the 
limit of zero momentum transfer (the 'standard' cross- 
section), F^{q) must be factored into ([T|) to describe the 
interaction with a particular nucleus. The definition of 
Ca in justifies the use of a normalised form factor 
F'^{q)=S{q)/S{0), where S{q) folds the spin structure 
functions Sjj,{q) with an arbitrary neutralino composi- 
tion in the isospin convention. The related functions 
Fppf (u) derived with Bonn CD were converted into the 
same Sjj,{q) normalisation used with Bonn A |il2i] by em- 
ploying Eq. (18) in Ref. [l3| and then parametrised with 
a 6*''-order polynomial for low q. As noted in Ref. the 
new calculation gives smaller spin structures than those 
found for simpler models: by a factor of ~2 for Xe^^^ 
and ~4 for Xe^'^^ at low q. More significantly, the spin 
factors are also smaller, in particular (Sp). 



TABLE I: Xe isotope parameters: nuclear spin J, isotopic 
fraction (nat. abundance), effective exposure and spin factors 
(Sn.p) with Bonn A 12] and Bonn CD 9] nn potentials. 







Bonn A 


Bonn CD 


J % (n.a.) 


kg-days 


{Sn} (Sp) 


(Sn) {S,} 


i29Xe 1/2 29.5 (26.4) 


37.7 


0.359 0.028 


0.273 -0.0019 


^=*^Xe 3/2 23.7 (21.2) 


33.7 


-0.227 -0.009 


-0.125 -0.00069 



Calculating the cross-section per nucleon allows com- 
parison of different target materials and with theoretical 
model predictions, which are usually computed for free 
protons and neutrons. This conversion is not straight- 
forward for the SD case since the cross-section has con- 
tributions from both proton and neutron terms, as in- 
dicated explicitly in ([2]). In addition, F^{q) is similarly 
contaminated by both couplings. A simple approach ex- 
ists that allows a straightforward calculation in a model- 
independent way by assuming that WIMPs couple pre- 
dominantly to protons or neutrons [3]. By setting, in 
turn, ttn—O and ap—0, no assumption is required as to 
the neutralino composition, either explicitly in the stan- 
dard cross-section or in the form factor. In this instance 
one may write, in the limit of zero momentum transfer: 



O-A 
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(3) 



where /J,p,„ is the WIMP-nucleon reduced mass. Once 
a A is obtained from experimental limits on the allowed 
nuclear recoil spectrum, ^ is used to calculate the corre- 
sponding SD cross-section for each isotope, and a}^^ . 
These are combined to obtain the total cross-sections: 

ZEPLIN-III is a two-phase xenon time projection 
chamber operating 1100 m underground at the Boulby 
mine. It is able to discriminate between nuclear recoil 
signals and the more prevalent electron recoil background 
by measuring both the scintillation light (SI) and the ion- 
isation charge (S2) produced by interactions in its 12 kg 
liquid xenon target. 3D reconstruction of the interaction 
vertex, a strong electric field in the liquid (3.9 kV/cm) 
and a geometry which avoids surfaces near a central fidu- 
cial region allow powerful discrimination down to low en- 
ergies. Data analysis procedures are detailed with our SI 
result 0] and briefly summarised here. The instrument 
is further described in Refs. [Hill]. 

A WIMP acceptance region was defined in the 
range [10.7,30.2] keV recoil energy and [/i-2cr,^] in the 
log]^Q(S2/Sl) discrimination parameter, where ^ and a 
describe the means and standard deviations of log- normal 
distributions fitted to the nuclear recoil population pro- 
duced by elastic neutron scattering over that energy 
range. An effective exposure of 128.7 kg-days was accu- 
mulated during the 83-day science run between 27**^ Feb 
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and 20* May 2008. This exposure is derived from a 
fiducial mass of 6.5 kg, defining a total 'geometrical' ex- 
posure of 454 kg-days, and a factor subsuming all energy- 
independent hardware and software inefficiencies (and 
the restricted acceptance). An energy-dependent detec- 
tion efficiency, which reaches unity near 14 keV recoil 
energy, is applied separately. Conversion between visible 
and nuclear recoil energies utilises the varying quenching 
factor (QF) discussed in Ref. 0. The exposure of the 
odd-neutron isotopes reflects their relative abundance. 
Our xenon was depleted from high-mass isotopes, es- 
pecially ^'^^Xe used in 0!//3/3-decay experiments. This 
enhances slightly the isotopic composition in ^^^Xe and 
^^^Xe relative to natural abundance (see Table I), as con- 
firmed by residual gas analysis using mass-spectroscopy. 

Seven events were observed in the acceptance region, 
all near the upper boundary of the discrimination pa- 
rameter. These were shown to be statistically consistent 
with the tail of the electron-recoil background popula- 
tion. To test this hypothesis in the presence of a hy- 
pothetical WIMP signal, a maximum likelihood analy- 
sis was performed on the acceptance region. This re- 
turned a most likely signal of zero WIMP events, with 
90% upper limits increasing from 2.45 to 3.0 events with 
increasing WIMP mass. The experimental limit on the 
WIMP-nucleus cross-section, a a, was calculated as de- 
scribed Ref. [l3|. We assumed an isothermal dark mat- 
ter halo with truncated Maxwellian velocity distribution 
with characteristic velocity tio=220 km/s, galaxy escape 
velocity Wesc^SOO km/s. Earth velocity W£;=230 km/s 
and WIMP mass density po=0.3 GeV/cm^. The limit 
on the differential recoil spectrum was corrected by the 
normalised nuclear form factor and detector parameters 
such as energy resolution and detection efficiency. 

The SD cross-section limits, with no assumption on the 
coupling strength to neutrons and protons, are shown 
in Figure 1 for the two spin structures. Values at 
55 GeV/c^ WIMP mass, the minimum of the curves, are 
given in Tabic II. Our Bonn A result surpasses that from 
XENONIO 18] above 100 GeV (taking into account new 
QF measurements for xenon pjj ) . Together these experi- 
ments place the world's most stringent limit on (t„, when 
this nn potential is assumed. However, Bonn CD affects 
these limits unfavourably: cr„ is just under twice higher, 
but Up increases by orders of magnitude, virtually elim- 
inating the sensitivity to WIMP-proton scattering (this 
curve is not shown in the figure as it would fit poorly 
within the range plotted). Naturally, these corrections 
apply to all xenon detectors. We note that comparison 
with other experiments is not always straightforward: 
i) different spin structures may be used for the same 
isotopes, as new calculations are still emerging; ii) the 
model-independent approach of Ref. [l3| is not used uni- 
versally (e.g. a combined SI and SD limit is extracted in 
HDMS, where a novel background subtraction technique 
is also employed [20| ) and iii) statistical significance may 
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FIG. 1: Upper limits on pure WIMP-neutron and WIMP- 
proton SD cross-sections. In addition to ZEPLIN-III with 
both nn potentials, we show other xenon experiments in 
black: ZEPLIN-I ZEPLIN-II [H and XENONIO (Bonn 
A) Additional curves are CDMS-IlJi|, COUPP [2I, 

EDELWEISS J2I, KIMS H^j, NAIAD [2||7pICASSO [H] 
and SIMPLE [30] . The pure-proton indirect limit from Super- 
Kamiokande is also shown [sj]. The DAMA evidence region 
interpreted as a nuclear recoil signal in a standard halo [Tj is 
indicated in green. The hatched area is the tip of the 95% 
probability region for neutralinos in the Constrained Minimal 
Supersymmetric Standard Model (CMSSM) |3^. 

vary, as with the XENONIO resuh, for which a CL<90% 
is ackowledged for the upper limit (c.f. note ]26] in ^21]). 

The allowed region of a„ — parameter space can be 
derived from the experimental cross-section limits jl3 ]: 

5:(^±f=)"<55^ (4) 

i VVCTp* Vcr„V ^^GtfA^p 

where i labels the two Xe isotopes, and the sign in paren- 
thesis is that of (Sp) / (Sp) . For xenon and other multi- 
isotope targets, Q defines an ellipse - albeit an elon- 
gated one - which reduces to two parallel lines for single- 
isotope experiments. The region allowed by each ex- 
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TABLE II: Spin-dependent cross-section limits (pb) at Mw^55 GeV for ^^^Xe and ^^^Xe, and the combined ZEPLIN-III result. 
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FIG. 2: Allowed regions in the a„ — ap parameter space for 
50 GeV WIMP mass. 



periment lies within the corresponding ellipse. Figure 2 
shows our Bonn CD result (nearly vertical lines, reflect- 
ing the poor constraint on Up) for a reference WIMP mass 
of 50 GeV/c^ (the cross-sections coincide with those in 
Table II within the precision shown) . The Bonn A result 
(not shown) is very similar to that from XENONIO. 

In conclusion, new experimental limits on the SD 
WIMP-nucleon cross-section are placed by the first sci- 
ence run of ZEPLIN-III operating at Boulby. A fiducial 
exposure of ^450 kg-days revealed a most likely signal of 
zero events, which allow 90% CL exclusion of the pure 
WIMP-neutron cross-section above (t„=1.8x10-^ pb for 
55 GeV/c^ WIMP mass. New spin structure calculations 
based on the Bonn CD nn potential were used for ^^^Xe 
and ^'^^Xe. These increase the cross-section limits rela- 
tive to previous calculations - quite dramatically in the 
case of the proton. Xenon targets are still the very sen- 
sitive to the WIMP-neutron interaction. Our result adds 
weight to the exclusion of the DAMA evidence region for 
a 50 GeV/c^ WIMP causing nuclear recoils. Theoret- 
ical predictions in the CMSSM [s^l point to x-nucleon 
cross-sections below ^SxlO""^ pb near 100 GeV mass at 
95% CL. ZEPLIN-III, XENONIO and CDMS-II provide 
the leading constraint on SD x-neutron scattering at that 
mass, nearly missing that probability boundary. Experi- 
mental efforts are still far from probing the favoured pa- 
rameter space for the x-proton interaction (the CMSSM 
best- fit regions suggest ~ 1 x lO-"* pb for both nucleons). 
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